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We present a complete calculation of the contributions to the effective leptonic weak mixing angle, 
sin 2 ^| ! g >t , generated by closed fermion loops at the two-loop level of the electroweak interactions. 
This quantity is the source of the most stringent bound on the mass, Mh, of the only undiscovered 
particle of the Standard Model, the Higgs boson. The size of the corrections with respect to known 
partial results varies between —4 x 10 -5 and —8 x 10~ J for a realistic range of Mh from 100 to 300 
GeV. This translates into a shift of the predicted (from sin 2 alone) central value of Mh by +19 
GeV, to be compared with the shift induced by a recent change in the measured top quark mass 
which amounts to +36 GeV. Our result, together with all other known corrections is given in the 
form of a precise fitting formula to be used in the global fit to the electroweak data. 
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The search for the Standard Model Higgs boson lies 
among the most important objectives of present elemen- 
tary particle physics. The experimental discovery will be 
possible at the Large Hadron Collider (LHC) within a 
mass range reaching up to 1 TeV. On the other hand, 
it is more than desirable to have as stringent indirect 
bounds on Mji as possible with the help of precision 
measurements. Should the Higgs boson be discovered, 
these bounds will serve as a strong test of the model. 

In this letter we study the quantity that has the high- 
est weight in the combined fit to electroweak data as far 
as Mjj prediction is concerned, which is the effective lep- 
tonic weak mixing angle, sin 2 ■ It can be defined 
through the form factors at the Z boson pole of the ver- 
tex coupling the Z to leptons (I). If this vertex is written 
as i h^igv - gAls)l Z» then 

sin 2 = 1/4(1- Re (g v /g A )). (1) 

At tree-level this amounts to the sine of the weak mix- 
ing angle sin 2 0w = 1 — M w /M^ m * ne on-shell scheme. 
In practice, sin 2 fl^g^is derived from various asymmetries 
measured around the Z boson peak at e + e~ colliders af- 
ter subtraction of QED effects. The current experimental 
value is 0.23150 ± 0.00016 Q. The high precision quoted 
and the expected size of the radiative corrections make 
the result indispensable for a precise prediction of Mh. 
A lot of effort has been put into the theoretical calcula- 
tion of sin 2 O^q 1 - Besides the one-loop contributions also 
higher-order QCD corrections 0, 0] are known. How- 
ever, for the electroweak two-loop corrections, only the 
leading term in the large Mh expansion and the lead- 
ing || and subleading Q terms in the large top quark 
mass expansion are available up to now. The goal of the 
present work is the calculation of the complete two-loop 



electroweak contributions with one or two closed fermion 
loops. 

The prediction Eq. (JTJ does not use My/ as input pa- 
rameter, but the results are given by using the very pre- 
cise measurement of the Fermi constant, G^, from the 
muon decay lifetime to derive My/ ■ Consequently the cal- 
culation of sin 2 as a function of G M involves also the 
computation of the radiative corrections to the relation 
between G M and Mw ■ For the electroweak two- loop cor- 
rections with closed fermion loops considered here, this 
has been carried out in Ref. Q. We will, therefore, also 
use the quantity Ak, 

sin 2 ^ c ff pt = (1 - M W /M 2 Z ) (1 + Ak) , (2) 

which is only weakly sensitive to My/, but encompasses 
the loop corrections to the Z form factors. 

We focus in this letter on the discussion of our main 
results. A detailed description of the calculation will be 
given in a forthcoming publication. Here, we note only 
that the contributions to the form factors can be divided 
into two major parts. The first one comprises the terms 
from renormalization. We use the on-shell renormaliza- 
tion scheme, similarly to the previous calculation of My/ 
yfl. The second one consists of the bare two- loop ver- 
tex diagrams, the total number of which approaches five 
hundred. Upon restriction to those containing a closed 
fermion loop we count only a few tens, which can be 
cast into four topologies as shown in Fig. ^ There is no 
dependence on the Higgs boson mass in the pure two- 
loop vertex diagrams, since CP conservation makes d) 
vanish. It is convenient to subdivide the remaining dia- 
grams into those containing a top quark line and those 
containing only light fermion lines. The former can be 
evaluated with the large top quark mass expansion, us- 
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FIG. 1: Two-loop vertex diagrams entering the calculation. 

ing the smallness of the ratio M§/m| « 1/4. We have 
convinced ourselves that an expansion up to (M^/m^) 5 
is sufficient to obtain an intrinsic precision of the order 
of 10~ 7 . The diagrams with only light fermion lines in- 
troduce also an important simplification: they have just 
two scales at most, Mw and Mz, since at the level we are 
considering, light fermion masses can be safely neglected. 
The problem thus contains only one variable, and lends 
itself naturally to the approach of differential equations 
||. A prerequisite for this method is the complete re- 
duction of the integrals to a small set of independent 
masters. This has been achieved with the C++ library 
DiaGen/IdSolver and been checked for a number of 
diagrams by an independent calculation. At the end we 
obtained analytic expressions for all of the integrals but 
one. The latter, corresponding to diagram b) of Fig. Q 
has been evaluated by a one dimensional integral repre- 
sentation. All integrals have been checked by different 
expansions in physical and unphysical regimes and by 
numerical integrations based on dispersion relations [lflj 
and Feynman parameterizations til) . 

An interesting problem connected to two-loop vertex 
diagrams is the treatment of the 75 matrix in triangle 
fermion loops. We used the naive dimensional regular- 
ization with a four-dimensional treatment of resulting ep- 
silon tensors as already explained in Q- We observed, 
however, that the contributions are divergent due to the 
soft-collinear behavior of the diagrams with external on- 
shell massless fermions. This would undermine the cor- 
rectness of the approach if the dimension of space-time 
were the only regulator. We decided to use a finite pho- 
ton mass as the regulator at the expense of a subse- 
quent difficult expansion corresponding to a mixed Su- 
dakov/threshold regime. The difference between the full 
result and the result which would be obtained if all traces 
containing a single 75 were set to zero will be denoted in 
what follows with a tr7s subscript. 

We shall now discuss the numerical effect of the new 



TABLE I: Input parameters with errors where relevant for 
the present analysis. 



input parameter 


value 


M w 


80.426 GeV 


M z 


91.1876 + 0.0021 GeV 


Tz 


2.4952 GeV 


m t 


178.0 ± 4.3 GeV 


m b 


4.85 GeV 


Aa{M%) 


0.05907 ± 0.00036 


a s (M z ) 


0.117 + 0.002 




1.16637 x 10 -5 GeV~ ; 



two-loop result for the effective weak mixing angle. We 
focus on the contributions to Ak, Eq. (J2J), taking the 
current experimental value for Mw as input. The as- 
sociated error is not relevant for the analysis, since the 
final prediction uses G M as input, combining the radia- 
tive corrections to Mw and Are. We use the parameter 
values given in Tab. [I] Note that the experimentally de- 
termined W and Z boson masses correspond to a Breit- 
Wigner parametrization with a running width and have 
to be translated to the pole mass scheme used in our cal- 
culation , resulting in a downward shift . For Mw 
and Mz, this shift amounts to about 27 and 34 MeV, 
respectively. 

Tab. |n] contains the values of the one- and two- loop 
electroweak corrections in comparison with different com- 
ponents with a single fermion loop for different values of 
the Higgs boson mass. The full two-loop result in the 
third column corresponds to the sum of the fourth, fifth 
and sixth column plus the contributions with two fermion 
loops as well as the effect induced by the running, Aa, 
of the fine structure constant. In the one-loop result we 
have kept a finite b quark mass, which has an impact of 
the order of —4.5 x 10 -5 . The perturbative expansion is 
performed in the fine structure constant, a, and not in 
G M , since we want to avoid any uncontrolled "resummed" 
terms. The first observation is that the third quark fam- 
ily contributions are very large, which is expected, since 
they include the leadin g top -bottom mass splitting effects 
in the p parameter, Ap [13j. We have convinced ourselves 
that the result has the correct behavior for large m t Hi- 
lt is interesting that even though the light fermion contri- 
butions in Tab.lrfldo not contain the running, Aa, of the 
fine structure constant, they are sizable. The last column 
gives the values of the tr75 contribution. It has to vanish 
for vanishing mass splittings in the fermion families and 
can be at most logarithmic for large top quark masses, 
which explains its smallness. For small Mh, the total 
two-loop result is rather small, but we note that this is 
due to a fragile cancellation strongly dependent on mt- 
With the older value mt = 174.3 GeV, the result would 
be of the order of 5 x 10~ 4 for all values of Mh- 
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TABLE II: One-loop and fermionic two-loop electroweak con- 
tributions to Ak with Mw as input parameter. The sub- 
scripts "tb", "If" and "tr75" correspond to the contributions 
of single loops of the third quark family, of the light fermions 
(without the running of the fine structure constant) and of 
the "tr75" effects in the triangle fermion subloops (see text). 



M H 


0(a) 


C0 2 )ferm 




C9(a 2 ) lf 


0(« 2 )tr 75 


[GeV] 


x 10~ 4 


X 10~ 4 


x 10~ 4 


x 10~ 4 


x 10~ 4 


100 


438.94 


-0.63 


-16.96 


-2.84 


0.27 


200 


419.60 


-2.16 


-17.10 


-3.08 


0.27 


600 


379.56 


-5.01 


-16.89 


-3.77 


0.27 


1000 


358.62 


-4.73 


-14.90 


-4.25 


0.27 



TABLE III: Various QCD corrections to Ak and the only 
known pure three-loop electroweak irreducible contribution, 
stemming from Ap, in comparison with three-loop reducible 
effects. The input parameter is Mw- 



M H 


0(aa 3 ) 


0(aa 2 s ) 


0(a 2 a s m$) 


0(a 3 m e t ) 


reducible 


[GeV] 


x 10" 4 


x 10" 4 


x KT 4 


x 10" 4 


x 10" 4 


100 


-36.83 


-7.32 


1.25 


0.17 


0.92 


200 


-36.83 


-7.32 


2.08 


0.09 


0.94 


600 


-36.83 


-7.32 


4.07 


0.07 


0.97 


1000 


-36.83 


-7.32 


5.01 


0.99 


0.98 



In order to provide the most precise prediction for 
sin 2 O^^m the SM we must use the muon decay con- 
stant, Gfj,, as input parameter. The procedure to derive 
Mw from G M is described in detail in ^4|. In analogy 
to that work, we do not want to perform any "resumma- 
tions". Instead, we include both in Ar and sin 2 fljrj* all 
known effects in expanded form. Besides the electroweak 
two-loop terms presented above, these effects encompass 
QCD corrections to the one-loop prediction at the two- 
Q and three-loop level and also the recently obtained 
0(a 2 a s ra\) and 0(a 3 m|) corrections to Ap 15]. We 
kept again a finite b quark mass in the 0{aa s ) correction, 
which has an impact of4.5xl0 -5 , almost completely can- 
celing the similar effect in the O(a) prediction. Consis- 
tency requires that we also take leading reducible effects 
at 0(a 2 a s ) and 0(ev 3 ) into account. It turns out that 
separate terms as e.g. Cy/ / Sy/ApAa 2 are quite sizable, 
but when summed cancel each other as seen in Tab. Mil 
We stress once more at this point that the same effects 
have been included in Ar and in sin 2 # 1< 3 >t . This means 
in particular that, contrary to |14| . we do not take the 
bosonic corrections |16j to Ar into account. Such precau- 
tions are enforced by the sensitivity of sin 2 0^ t to My/, 
since a 1 MeV shift in the latter causes a shift of about 
—2 x 1CP 5 in the former. 

Our complete result is summarized by the following 
fitting formula, which reproduces the exact calculation 



with maximal and average deviations of 4.5 x 10 -6 and 
1.2 x 10~ 6 , respectively, as long as the input parameters 
stay within their 2a ranges and the Higgs boson mass in 
the range 10 GeV < M H < 1 TeV, 



sm^ff 



d 2 L 2 



s + d\Lu -r u, 2 u H 
+d 5 A a + d 6 A t + d 7 A 
+d 9 A as + d w A z , 



d 3 L% + d 4 (A 2 H - 1) 



■d 8 Ai(A H -l) 



with 



L H = log 



M H 



100 GeV 



A 



M, 



H 



H 



Aa 



0.05907 
a s (M z ) 
0.117 



-1, 



A z = 



100 GeV 

m t 



(3) 



(4) 



178.0 GeV 

M z 
91.1876 GeV 



- 1, 



and 



s = 0.2312527, 
di = 4.729 x 10~ 4 , 
d 3 = 3.85 x 10~ 6 , 
d 5 = 0.0207, 
d 7 = 1.82 x 10~ 4 , 
d 9 = 3.98 x 10~ 4 , 



d 2 = 2.07 x 10~ 5 . 
di = -1.85 x 10" 
d 6 = -0.002851, 
rf 8 = -9.74 x 10" 
d m = -0.655. 



(5) 



The impact of this result is shown in Tab. IIVI where 
we compare our prediction with the previous result as 
given in the fitting formula in |17| and implemented in 
ZFITTER [13]. The difference varies from roughly -4 x 
10~ 5 to -8 x 10~ 5 for the M H range from 100 to 300 
GeV, which is the preferred mass region inferred from 
precision electroweak data. These values reach half of the 
experimental error and induce an important shift in the 
central value of Mh derived from sin 2 alone. With 
the most recent value of the top quark mass given in 
Tab.njthe result shifts the central value from 149 GeV to 
168 GeV, to be compared with the shift induced by the 
new mt measurement which gives a jump from 132 GeV 
to 168 GeV. The formula Eq. JHJl has been implemented 
in the most recent version of ZFITTER, version 6.40 [lfSj . 

Besides providing an up to date fitting formula, it is 
necessary to discuss the error on the theoretical predic- 
tion connected with the unknown higher order contri- 
butions. Here one has to incorporate the treatment of 
the error of the My/ prediction, since the final prediction 
for sin 2 ^^ )t takes G M as input. In particular, there are 
some cancellations between the radiative corrections to 
My/ and the Z decay form factors that go into sin 2 . 
We take the point of view that these cancellations are 
natural and discuss both quantities in conjunction. To 
this end, we use geometric progression from lower orders 
to estimate the missing higher-order contributions and 
add them quadratically at the end. In units of 10~ 5 we 
assign the following errors: corrections of 0(a 2 a s ) be- 
yond m\ vary between 2.3 and 2.0 for Mjj between 10 
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TABLE IV: Difference between the result Eq. and the 
previous result including terms of 0(a 2 m 2 ) from |g], obtained 
from the ZFITTER implementation (left column) or from the 
fitting formula from [17fl. 



M H 
GeV 


(A sin 2 r? lcpt ) 

V / ZFITTER 

xhT 4 


xlO" 4 


100 


-0.45 


-0.40 


200 


-0.69 


-0.72 


300 


-0.85 


-0.83 


600 


-1.17 


-0.94 


1000 


-1.60 


-1.28 



0.2325 



0.232 



0.2315 



exp. lower bound MH = 1 14.4 GeV 




300 500 700 

M„ [GeV] 



FIG. 2: The sin 2 6 e £^ prediction against the current experi- 
mental value, with la bands from the experimental input. 



GeV and 1 TeV, corrections of 0(a 3 ) between 1.8 and 
2.5, 0(aal) between 1.1 and 1.0, 0(a 2 a 2 s ) between 1.7 
and 2.4, and finally the bosonic corrections at 0(a 2 ) are 
expected to be of the order of 1.2. This gives an error 
varying between 3.3 and 3.5, and we take as our esti- 
mate the latter, largest, value. To account for possible 
deviations from the geometric progression of the pertur- 
bation series, we included an additional overall factor of 
a/2, giving a final error of 4.9 x 1CP 5 . 

If we take into account all of the input parameter er- 
rors, our prediction can be compared with the experi- 
mental value as shown in Fig. [3 

In conclusion, we have calculated the complete 
fermionic corrections to sin 2 6 ) ^ t at the two-loop level 
and obtained a sizable contribution when compared to 
the previously known leading and subleading terms in 
the top quark mass expansion. Together with our result 
for the W boson mass and recently obtained three-loop 
terms, we are able to give the most up to date prediction 
to be used in the global fit to electroweak data. Further- 
more, we implemented the result in the program ZFIT- 
TER, widely used for this purpose. 
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